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ABSTRACT 

We explore the effect of the magnetic field when using realistic three- 
dimensional convection experiments to determine solar element abundances. By 
carrying out magnetoconvection simulations with a radiation-hydro code (the 
Copenhagen stagger code) and through a-posteriori spectral synthesis of three 
Fe I lines, we obtain evidence that moderate amounts of mean magnetic flux cause 
a noticeable change in the derived equivalent widths compared with those for a 
non-magnetic case. The corresponding Fe abundance correction for a mean flux 
density of 200 G reaches up to ~ 0.1 dex. These results are based on space- and 
time-averaged line profiles over a time span of 2.5 solar hours in the statistically 
stationary regime of the convection. The main factors causing the change in 
equivalent widths, namely the Zeeman broadening and the modification of the 
temperature stratification, act in different amounts and, for the iron lines con- 
sidered here, in opposite directions; yet, the resulting |Aloge0(Fe)| coincides 
within a factor two in all of them, even though the sign of the total abundance 
correction is different for the visible and infrared lines. We conclude that mag- 
netic effects should be taken into account when discussing precise values of the 
solar and stellar abundances and that an extended study is warranted. 
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abundances — Sun: granulation — Sun: photosphere 



Introduction 



The advent of increasingly r ealistic three-dimens i onal convection s imulations based on 



radiation-hydro dyn amics codes (jPravins et al.l Il98ll : iNordlundl Il985l : iFreytag et al.l 12002 



Vogler et al.l |2005| ) has led to the conclusion that horizontal inhomogeneities and depar- 
tures from Local Thermodynamic Equi librium ( LTE] are crucial for an accurate abun- 
dance determination (see the review by lAsplundl 120051 ). Three-dimensional hydrodynamic 
(HD) models have allowed to remove many uncertainties in abundance determination, in 
particular those comi ng from the use of micro- and macroturbulence velocity parameters 
(lAsplund et al.ll2000al ). Yet, several sources of uncertainty remain, resulting, e.g., from pos- 
sible blending of lines, errors in continuum measurements, equivalent width determination 



(IHolweger. 



^ock. k Bardlll995l: 



B 



ackwell. Smith, fc Lynas-Graylll995l : lKostik. Shchukina fc Rutten 



1996| : lAvresl l2008aljbl:l Caffau et al.lboosi ) or calculation of collisions with H atoms (jPabbian et al 



20061 . l2009l ). Another potentially important source of uncertainty stems from the presence 



of significant amounts of magnetic flux in the photosphere of the quiet Sun 



Sanchez Almeida. Emonet fc CattaneolboOS : Shchukina fc Trujillo Bueno 



2003 



Stenfio 



1982 



20041) and its possible influence on the average width of the line proflles (jStenflo fc LindegrenI 



Truiillo Bueno. Shchukina 



1977i ). The purpose of this Letter is to use 3D magnetoconvection (MHD) models to esti- 
mate the influence of the magnetic fleld on the spectral proflles of a number of Pel lines, 
and to translate the difference between the equivalent widths in MHD and HD models to 
the corresponding Fe abundance corrections. 

The impact of magnetic flelds on chemical abundance determinations has so far been 
widely neglected. It was implicitly assumed that the magnetic fleld should not play a ma- 
jor role in changing the sh apes of spatially-averaged spectral lines, as these averages are 
weighted towards granules ( Nordlund. Stein. &: Asplund 20091), while magnetic con centra- 
tions are thou ght to reside mainly in intergranular lanes ( jKhomenko et al" 2003). Still, 
Borrerd (120081 ) found that the neglect of magnetic broadening can lead to a systematic error 
of up to 0.1 dex in the abundance determination of iron. The author used ID spectral synthe- 
sis with uniform magnetic fleld in a plane-parallel atmosphere, therefore inevitably neglecting 
the effect of the fleld on the temper ature stratiflcation. However, t he latter effect, although 
indirect, can be important (see, e.g. Socas-Navarro fc Norton 2007|) and must be considere d 
alongside the direct effect through Zeeman broadening (see, e.g. IStenflo fc LindegrenI 119771 ). 



Iron is one of the crucial chemical elements for a variety of reasons like: (a) it is widely 
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used to define a metallicity scale, hence as a proxy for tlie overall metal content; (b) tlie effect 
of tlie magnetic field on the line pr ofiles seem to be more pronounced for this element than 
for, e.g. Si, C and O ( lBorrerdl2008l ): (c) the large amount of iron lines in the solar spectrum 
allows to consider features with high Lande factor. It is therefore advisable to start with 
iron when studying the magnetic effects on abundance determinations from MHD models. 
In this Letter, we present results based on a series of three-dimensional magneto-convection 
simulations. Our aim is to explore, understand and constrain the possible impact on the 
solar abundance determination of self-consistent, simultaneous 3D and magnetic effects. 
These effects are studied via comparison of the equivalent widths for the magnetic cases 
with those obtained in HD conditions. To derive the necessary abundance correction, we 
match the MHD equivalent width with that obtained on the basis of HD runs with changed 
abundance. Our analysis focuses on quiet Sun conditions. We conclude that magnetic effects, 
in particular via modification of the average temperature stratification, should be taken into 
account when considering the solar abundance problem. 



2. Simulations and spectral synthesis 



The resul ts in this paper are based on s i mulations we carried o ut using the Copenhagen 



stagger-code (iNordlund fc GalsgaardI 119951 : ICarlsson et al.l |2004J : IStein &: Nordlundl 12006 



De Pontieu et al.l |2006[ ). a high-order finite-difference Cartesian MHD code using hyper- 



diffusivities. All calculations reported here are for a domain of 6 x 3 x 6 Mm^, with a grid of 
252 X 126 X 252 points and uniform spacing in the horizontal directions. The vertical direc- 
tion is non-uniform and reaches a maximum resolution of ~15 km in the pho tosphere. This 



grid sh ould provide sufficient resolution for the purpose of abundance studies (lAsplund et al. 
2000bh . 



To initiate the MHD runs, a vertical uniform magnetic field of strength Bo was intro- 
duced into an already evolved HD snapshot. An HD run and three different MHD cases 
(with Bq = 50, 100, 125 G) were evolved for up to several solar hours, with snapshots taken 
every 30 seconds. When the implanted magnetic field is sufficiently strong, as the simula- 
tion evolves the continuum intensity starts to show fine bright threads /points and fragments 



(Fig. [H), as in high - resolution solar observations (e.g. ISanchez Almeida et al.l l2004l 12010 
Shelvag et al~ll2004l ). 



To simplify the task, the spectral synthesis was performed in LTE. We considered the 
three lines of Fel at A = 608.27, 624.07 and 1564.85 nm, which have EPL = 2.223, 2.223 
and 5.43 eV; log gf = -3 572, -3.390, -0-670 and ges = 2.00,0.99 and 2.98, respectively. 
We used the LILIA code ( ISocas-Navarrd 120011 ) . and set the resolution to 1.28 pm over 200 
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wavelength points. 

The synthetic spectra were computed using no enhancement factor to the Van der Waals 
broadening formula in the spectral synthesis. To obtain the final synthetic profiles, emergent 
spectra were calculated for each vertical column in the grid. Averages were then calculated 
for each wavelength over the horizontal directions and in time, using 95 solar minutes of the 
stationary stage of the simulations. 



3. Results 



3.1. Goodness of the simulations 



In Fig. [21 we show (solid lines) the spatially- and temporally-averaged synthetic line 
profiles from our HD and 200 G experiments (top panels), and the differences in intensity 
between the spectral profiles derived in the MHD cases and those obtained for the HD case 
(bottom panels). Each of the profiles shown was normalized to its own continuum level. 



Observational data taken from the FTS spectral atlas (iBrault et al.lll987l ) are also shown 
in the figure for visual comparison (top panels, empty circles). The fit of the line profiles 
is satisfactory, with the MHD models performing better than the HD one. In any case, 
since the objective is a differential analysis between MHD and HD, the scope of the work 
presented here does not require obtaining extremely high accuracy in the fitting of observed 
line profiles. Remaining discrepancies in the fitting, especially in the line cores, can be 
explained by the higher temperatures in the relevant layers of our MHD models, which will 
lead to stronger neutral iron overionisation and are thus likely to induce larger departures 
from the LTE assumption adopted here. 

More importantly with regards to our aim here, we verified that our equivalent width 
calculation converged with regards to the number of sampled spectral wavelengths. More- 



found in the literature (IGurtovenko & KostikI 


1986; 


Blackwell. Lvnas-Grav. & Smith 


1995 


Holweeer. Kock. &; BardI 


I995I; 


Kostik. Shchukina &; Rutten 


1996). 



A further check of the goodness of our simulations is possible: we obtain an excellent 
match be tween our continuum intensity val ues and those expected from the literature (see 
figure 2 of lTrujillo Bueno fc Shchukinall2009l ). For example, in units of erg/ (cm^ s ster A), we 
get Stokes I values of ~ 3.13 ■ 10^ (continuum around the Bell 313.04 and 313.11 nm lines), 
~ 4.88 • 10^ (continuum around the Fel 401.02 nm line), ~ 2.91 ■ 10^ (around the Fel 630.15 
and 630.25 nm lines) and ~ 2.43- 10^ (continuum around the Fel 700.06 nm line). The latter 
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two va lues match almost perfectly with the observational data of lAyres. Plymate. &: Keller 
( 120061 ) ■ On the other hand, our UV values, while still agreeing to a good level, may be 
slightly too high due to our neglect of line haze on continuum at those wavelengths. 

We have also used Stokes V polarization signature results (Fig. [H right panel) to derive 
Stokes V amplitudes. Depending on the wavelength considered, we find Stokes V average 
amplitudes (measured in units of Stokes I continuum intensity) of the order of 0.01 — 0.03% 
("100 G" case). Thi s corresponds, when calibrating our Stoke s V signals using the weak- field 



approximation, (e.g lMartfnez Gonzalez fc Bellot Rubidl2009t ). to an average unsigned verti 



cal magnetic flux of the order of 10 — 50 G. These numbers are typical for solar internetwork 
and network regions, as deduced when considering the effect of limited spatial resolution 
of observations. This result confirms that our simulations are representative of solar quiet 
regions, a description covering 90 percent of the solar surface. 



Regarding limb-darkening, recent work by lPereira. Asplund fc Kiselman I (j2009al ) shows 
that related uncertainties in the simulations play only a minimal role for solar abundances. 
Our differential approach ensures that the small impact of any discrepancy found against 
limb darkening measurements can be neglected here. 



3.2. Thermal versus magnetic effects 



Figure [3] compares the average temperature stratification of the different magnetic runs 
versus optical depth: the averaging was done on a fixed log(T5oo) scale, and over all snap- 
shots. As seen in the figure, in the relevant layers of formation of the absorption features 
we consider, (T) increases with increasing magnetic flux in the simulation s, a result ex - 
pected from radiation transf er considerations in strong solar magnetic elements ISpruitI (Il976l ): 
Schiiessler &: Solankil (119881 ). 



The formation heights of the relevant lin es, for the HP case, are: -0 .276 Mm for Fel 
624.07 nm and -0. 229 Mm for Fel 608.27 nm dgurtovenko fc Kostikl ll986h. and -0.135 Mm 
for Fe 1 1564.85 nm lShchukina fc Trujillo Buend (120011 ). which correspond to an optical depth 
of between ~ —1.6 and ~ —0.7. As seen in the figure, at these optical depths there is a 
significant temperature increase in the MHD cases compared to HD. The maximum temper- 
ature difference, arises at log(r5oo) ~ —2.5, reaching up to 125 K for the run with average 
vertical magnetic flux density of 200 G. Together with the direct effect due to the presence of 
a magnetic field, this has direct implications for the equivalent width of the average profiles 
and, thus, for the abundance corrections. 



As seen in Fig. [21 the line core tends to get weaker in the MHD cases for all of the Fe 
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lines considered, due to the increased average temperature of the corresponding atmospheric 
models (Fig. |3]). For the two visible lines, this indirect influence on the line formation turns 
out to be the dominant effect, with the direct effect of the magnetic field on the lines being 
small. Due to its IR wavelength and large magnetic-sensitivity, the 1564 nm line experiences 
instead significant direct Zeeman broadening, which, in terms of resulting equivalent width, 
goes in the opposite direction of the indirect (core-weakening) effect^ 



3.3. Differential abundance corrections 

To obtain a differential abundance correction between the HD and MHD cases, we 
modify the Fe abundance upward or downward of the standard value loge0(Fe) = 7.50 



( IShchukina fc Trujillo Buend 1200 ll : lAsplund et al.l l2009l ) in the calculation of the average 



spectra for the HD series. This is done in steps of 0.01 dex, until the corresponding theoretical 
equivalent width matches the one derived using the standard iron abundance when computing 
the average spectra for the MHD case of interest. This method allows us to translate the 
effects of the magnetic field into abundance corrections. 

The equivalent widths of the three spectral lines calculated with the standard Fe abun- 
dance are presented in Fig. H] (filled circles). The figure also gives (dashed horizontal lines) 
the equivalent width values for the HD case calculated with the iron abundance changed 
in steps of 0.01 dex. Since for each case we normalise the average profiles to the relevant 
continuum level, changes in equivalent width are due to actual changes in the shape and 
depth of the absorption features. By comparing the filled circles and the horizontal lines in 
each of the three panels of Fig. HJ we see that the maximum abundance difference is reached 
in all cases for the runs with the highest field strength, B = 200 G. The actual numerical 
values for the maximum difference are: |A loge0(Fe)| ~ 0.06 dex for Fel 1564.85 nm; ~ 0.09 
dex for Fel 624.07 nm and slightly above ~ 0.07 dex for Fel 608.27 nm. 

The variations of the equivalent width in Fig. H] are caused both by direct magnetic 
effects [i.e., via Zeeman broadening) and, indirectly, via changes in the stratification leading 
to different amount of absorption in the line core. As the name implies, Zeeman broadening 
always causes an increase in the equivalent width. In contrast, the indirect contribution 
resulting from the increased average temperatures in the MHD models is a hne-w eakening 



effect for the three chosen lines. This is to be expected (see discussion in lGraylll992l Chapter 



13) for neutral Fe lines in the solar photosphere, where iron is mostly ionized in the relevant 



^Note that this hne is formed significantly deeper, where the temperature increase compared to the HD 
case is less significant. 
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layers. One can estimate the relative importance of the direct and indirect effects by ad-hoc 
setting Bq = while keeping all other variables unmodified when doing the spectral synthesis 
for the MHD snapshots. The open circles in Fig. H] show the equivalent widths resulting 
from carrying out this experiment for all three lines. For the visible lines the indirect effects 
dominate, which can be expected since they form at heights in the photosphere where the 
temperature differences between the MHD and HD models become more important. The 
direct effect for those lines leads to a roughly linear increase of equivalent width with B, 
which is compatible with the Zeeman broadening origin of this effect. For the infrared Fel 
1564.85 nm line, in turn, the direct magnetic effects become dominant to the point that they 
reverse the sign of the abundance correction in that case. The sensitivity of this spectral line 
to the magnetic field is in fact very high, given its longer wavelength and the fact that its 
Lande factor is gi ~ 3. Also, for this infrared Fel line there seems to be a saturation effect, 
although no asymptotic horizontal regime is reached in the range of B used in the present 
study. 



Conclusions 



Pereira. Kiselman &: Asplund I (|2009[ ). iPereira. Asplund fc Kiselman I ( l2009bl ). ICaffau et al. 



Caffau et al. 


( 


2008 


), 


f2009b). 


Caffau et al. 



the agreement with solar-neighbourhood and local interstellar-medium abundances, the pro- 
posed large reduction of the solar overall metallicity has led to controversies, especially from 
the helioseismology community, because it spoils the previous agreement between interior 
sound speed profiles derived from helioseismic measurements and predictions based on the 
classical solar model teahcall et al. 2005; Antia &: Basu 2006 ). 



To settle the open issues, a thorough exploration of the effects on line formation of 
the interaction between magnetic field and convection is a necessary step. In this letter we 
have presented a study of that kind using three representative Fe lines. We have deter- 
mined that magnetic fields have a noticeable influence on the formation of iron lines and 
consequently on the solar Fe abundance determination. For average vertical flux densi- 
ties of 200 G, we obtain abundance corrections (to be applied in order to obtain the same 
equivalent widths in HD and MHD) between 0.06 and 0.09 dex in absolute value, while 
|Alogeo(Fe)| ~ 0.05 dex for all lines for the 100 G case. A flux density of order 10^ G, 
as in the range shown in Fig. |U is probably the relevant reference value since it has been 
shown to correspond to the order of the mean field strength of the photospheric plasma 
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( iTrujillo Bueno. Shchukina fc Asensio Ramosll2004[ ). 



Our results show that the direct (Zeeman) and indirect (stratification-related) effects of 
the magnetic field contribute in different amounts to the changes in the different lines. The 
indirect effects cause most of the equivalent width change of the visible Fe lines, with a small 
contribution from the Zeeman effect; for the infrared line, in contrast, the Zeeman effect 
is dominant, and in fact changes the sign of the abundance correction. The absolute value 
of the total correction is similar for all lines, so we can take it as the reference abundance 
change related with the presence or absence of magnetic field in the convection model. In 
fact, for the visible lines, the use of MHD models to match the observational data will lead to 
higher derived iron abundance than when using an HD model, since the predicted theoretical 
equivalent widths are smaller for the former cases than for the latter. The opposite applies 
for the Fel 1564.85 nm, since the increased temperature of the MHD models is more than 
compensated by the very large direct Zeeman line splitting. 

There are different avenues to explore with a view to confirming and extending all 
these results. O ne should consider further r eliable Fe I line pairs with different mag- 
netic sensitivity ( IVasilyeva fc Shchukinal I2OIOI ). and carry out an extensive study includ- 
i ng all Stokes parameters for the s ame spectral feature s as in previous 3D investigations 
( jShchukina fc Trujillo Buend I2OOII : lAsplund et al.l |2009| ) , but now based on MHD simula- 
tions. 



On the other hand, spectral synthesis using 3D convection m odels should aim to match 

key observational constraints such as spatial ly resolved line profiles ( jPereira. Kiselman fc Asplund 
2OO9I : iPereira. Asplund fc Kiselman Il2009bl ). center-to-limb continuum behaviour, calibrated 
continuum intensities at disk center, spectral energy distribution, and H lines. Finally, it is 
urgent to explore the magnetic effects in the Sun for lines of other crucial elements. 
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Fig. 1. — Emerging continuum intensity at 608 nm in the HD simulation run (left), and 
in the MHD run with "100 G" setup (center). Note the appearance of intergranular bright 
points after introducing magnetic flux. The right panel shows the amplitude of Stokes V, 
defined as max(|V|) in units of the continuum intensity, sign included. 
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Fig. 2. — Top panels: normalized theoretical profiles after spatial- and temporal-average 
for the Fel 608.27 nm (left), 624.07 nm (center) and (right) 1564.85 nm spectral lines, as 
derived from the synthesis of relevant spectral regions from the HD and MHD simulations. 
For the sake of clarity, MHD results are shown only for the simulation with highest magnetic 
flux. Spectral observational data from the Kitt Peak FTS Spectral Atlas are given for visual 
comparison (empty circles). Bottom panels: intensity difference between average normalized 
synthetic line profiles in HD and MHD. 
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Fig. 3. — Differences with the HD case in the temperature stratification (T) (with averaging 
taken at fixed optical depth, and over time), for i?o = 50 G (blue line), 100 G (green hne), 
and 200 G (red line) cases, plotted against a logio{T5oo) scale. The approximate photospheric 
region of formation for the lines studied is shown with a gray vertical band. 
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Fig. 4. — Filled circles: derived equivalent widths in the HD run and the three MHD 
runs, with separate panels for each spectral line. Empty circles: equivalent width in the 
experiments, when doing the spectral synthesis while setting B — (leaving the thermody- 
namics unchanged). The horizontal lines mark the equivalent width for the HD case using 
the standard solar Fe abundance (solid line), and for iron abundance changes in steps of 
|Aloge0(Fe)| = 0.01 dex. 



